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ABSTRACT 

We report observations of 15 high redshift (z = 1 — 5) galaxies at 350//m using the Caltech Sub- 
millimeter Observatory and SHARC-II array detector. Emission was detected from eight galaxies, 
for which far-infrared luminosities, star formation rates, total dust masses, and minimum source size 
estimates are derived. These galaxies have star formation rates and star formation efficiencies com- 
parable to other high redshift molecular emission line galaxies. The results are used to test the idea 
that star formation in these galaxies occurs in a large number of basic units, the units being similar to 
star-forming clumps in the Milky Way. The luminosity of these extreme galaxies can be reproduced 
in a simple model with (0.9-30) x 10 6 dense clumps, each with a luminosity of 5 x 10 5 L Q , the mean 
value for such clumps in the Milky Way. Radiative transfer models of such clumps can provide reason- 
able matches to the overall SEDs of the galaxies. They indicate that the individual clumps are quite 
opaque in the far-infrared. Luminosity to mass ratios vary over two orders of magnitude, correlating 
strongly with the dust temperature derived from simple fits to the SED. The gas masses derived from 
the dust modeling are in remarkable agreement with those from CO luminosities, suggesting that the 
assumptions going into both calculations are reasonable. 

Subject headings: galaxies: high-rcdshift — galaxies: starburst — infrared: galaxies — galaxies: ISM 
— galaxies: formation 



1. INTRODUCTION 

Studies of star formation in the early Universe are im- 
portant to an understanding of galaxy formation and 
evolution. Observations of quasar host galaxies and 
sub-millimeter galaxies with detectable molecular line 
emission offer an opportunity to extend such studies to 
high redshift (z > 2), albeit wi th a limited sample (see 
iSolomon fc Vanden BouTI 2005). Following that review, 
we refer to these galaxies collectively as EMGs (Early 
Universe Molecular Emission Line Galaxies). Several 
kinds of models of these objects have been proposed. 
Dust in quasar host galaxies could be heated by the ra- 
diation from the cent ral AGN (e.g. JGranato et ai1ll996l . 
lAndreani et al.l fl999f ) . While it is hard to rule out such 
models, general arguments t end to favor starb ursts as 
the main power source (e.g., IBlain e t al. 2002j). In ad- 
dition, detailed studies of well-known sources have in- 
dicated that star for mation dominates the contribution 
from the black hole dRowan- Robinson! [20001 IWeifi et al.l 
[2001 [Chapman et al.H2004 IWeifi et al.ll2007l ). 

There are also variations among models based o n star 
formation. lEfstathiou fc Rowan- Robinson! (|2003f) pro- 
posed a model in which the far-infrared emission arises 
from cirrus (relatively diffuse dust) heated by ultravi- 
olet photons leaking out from star formation regions. 
However, most recent analysts have focused on models 
in which the far-infrared radiation arises from dust that 
is intimately associated with a burst of star formation 
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(e.g., iNaravanan et al.ll2009aft . In the embedded star- 
burst models, the total luminosity of the dust continuum 
emission is a measure of the star formation rate, and 
the luminosity of the molecular line emission or the dust 
emission at long rest wavelengths measures the amount 
of material available for star formation. 

Here we report observations of the dust continuum 
emission at 350/im wavelength, obtained at the Cal- 
tech Submillimeter Observatory 4 . This observed wave- 
length falls roughly near the peak of the emission in 
the rest frame of the objects observed and is, there- 
fore, a desirable wavelength for observations to de- 
termine far-infrared luminosities. It will also provide 
stronger constraints on the characteristic temperature 
of the far-infrared radiation, which can distinguish be- 
tween cirrus models and models of embedded star for- 
mation, as the former predict cool (Tp < 30 K) dust 
(Efst athiou fc Rowan- Robinsod [20031 ). 
This work extends previous studies at 350Atm of high 



redshift molecular line emission galaxies ( Bcnfor d et al 
199i IWeifi et alJl2003t IBeelen et al. 11200a IWang et~aT 
2008al: I Wang et al.l I2008H ). Throughout this paper we 



have assumed a ACDM cosmology with Hp = 70 km s 
Mpc-\ fl m = 0.27, and fl A = 0.73 ([Spergel et al.|[2007T) . 

2. OBSERVATIONS 

Table 1 lists the objects we observed. They are, with 
two exceptions, objects previously unobserved at 350/im, 
chosen for the strength of their CO emission and the 
strength of their long- wavelength dust continuum emis- 
sion to maximize the success of detection. Table 1 also 
gives the coordinates observed, source redshifts from CO 
observations, the dates of the observations, the averaged 

4 The Caltech Submillimeter Observatory is supported by the 
NSF. 
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zenith angle during observation with atmospheric opaci- 
ties and integration times. 

Observations were conducted in several runs during 
2003 through 2007, using the Submillimeter High Angu- 
lar Resolution Camera II (SHARC-II) at the 10.4 m tele- 
scope of the Caltech Su bmillimeter Observa tory (CSO) 
at Mauna Kea, Hawaii (|Dowell et al J 120031) . SHARC- 
II is a background-limited camera utilizing a "CCD- 
style" bolometer array with 12 x 32 pixels. At 350 
/Ltm, the beam size is 8'.'5, with a 2'59 x 0.'97 field of 
view. Since the atmospheric transmission is very sensi- 
tive to the weather at the higher frequencies at which 
SHARC-II operates, most of our integrations were made 
at small zenith angles and under the best weather con- 
ditions at Mauna Kea, during which the opacity at 225 
GHz (t225GHz) was less than 0.06 in the zenith, which 
corresponds to an opacity of 1.5 at 350 ^m (see Table 1). 
For most of our obse rvations , the Dish Surface Optimiza- 
tion System [DSOS, ILeong I (l2005l) ] was used to correct 
the dish surface figure for imperfections and gravitational 
deformations as the dish moves in elevation. 

All the raw data were reduced with the Comprehen- 
sive Red uction Utility fo r SHARC-II, CR USH - version 
1.40a9-2 (|Kovacsll200"6al IKovacs l l2006bh We used the 
sweep mode of SHARC-II to observe all our sources. In 
this mode the telescope moves in a Lissajous pattern that 
keeps the central regions of the maps fully sampled. It 
works best for sources whose sizes are less than or com- 
parable to the size of the array, but causes the edges to 
be much noisier than the central regions, and can often 
result in noise at the edges that looks like real emission. 
To compensate for this, we used "imagetool," part of the 
CRUSH package, to eliminate the regions of each map 
that had a total exposure time less than 25% of the max- 
imum. This eliminates most, but not all, of the spurious 
edge emission. 

Pointing was checked every 1-2 hours during each run. 
The primary pointing sources were planets such as Mars, 
Uranus, and Neptune, and their moons, for example, Cal- 
listo. If no planets were available, we used secondary ob- 
jects such as CRL618 and IRC+10216. After averaging 
over all the runs, the blind pointing uncertainty is 2'.'1 
for azimuth and 3'.'1 for zenith angle. We corrected the 
pointing after each check, so these uncertainties actually 
represent upper limits. When reducing raw data with 
CRUSH, we also applied a pointing correction based on 
the statistics of all available pointing data (to remove the 
static error) and several pointings before and after ob- 
serving scans (to remove the dynamic error). This tech- 
nique improves the calculated flux densities of point-like 
sources and further reduces the pointing uncertainty 5 . 

To obtain the total flux densities of the sources in units 
of Jy, we have used Starlink's "stats" package to mea- 
sure the signal from targets in a 20" aperture (an aper- 
ture large enough to include all the emission from high-z 
galaxies) , and we measured the signal from calibrators in 
the same aperture. We used planets as calibrators when- 
ever possible, but some secondary calibrators were used 
when planets were not available. The Flux Conversion 
Factors (FCF) for an 20" aperture (C20) is defined to 
be the total flux density of a calibrator source in Jy di- 
vided by the signal in that aperture from the calibrator 

5 http:/ /www. cso. caltech.edu/~sharc/ 



in instrument units. Since CRUSH already includes an 
atmospheric correction, based on a fit to all calibrator 
observations during the night, the flux density within 
the 20" aperture is then obtained by simply multiply- 
ing the measured signal in the instrument units by C2o- 
The statistics of C20 for all the calibrators over all our 
runs indicates a calibration uncertainty of 20%. Detailed 
studies have been carried out of calibration uncertainties 
including uncertainties in the flux density of the calibra- 
tor, the airmass of the source, and differences in atmo- 
spheric opacity at different times throughout the night. 
These studies were based on surveys towards low-mass 
Galactic cores observed on the same runs as the galaxy 
observations. The results (Wu et al. 2007, Dunham et 
al. in prep), are that the 20% calibration error domi- 
nates over the airmass and atmospheric opacity errors. 
Therefore we take the systematic error to be 20%, which 
we add in quadrature to the random noise uncertainties. 

3. RESULTS 

Of the 15 objects observed, we have four clear detec- 
tions, (> 5<t), and four tentative detections (2a < S v < 
5a). These are listed in Table 2 with the detected flux 
density and noise in mJy, signal-to-noise ratio in a, and 
detection classification. The 350^tm images of the eight 
detected sources are given in Figures [T] and [21 where the 
circles indicate the region of reliable data reduction. 

We have collected data from the literature, focusing 
on long- wavelength photometric data. These data will be 
used to model the sources, both with simple graybody fits 
and with radiative transport models (3H|). The literature 
data are given in Table 3 and shown in Figure [3] 

Table 4 gives properties derived from the fits to the 
SEDs, corrected for magnification (noted as intrinsic 
properties) and redshift. We fit optically thin graybody 
spectra to the 350-1200/mi data points for each source, 
using spectral energy distributions (SEDs) of the form: 

S{v r est) = (1 + Z^SiVobs) 

= So ■ ( ^ Y +f> - (1) 

\\THz) exp(his res t/kTi)) — 1 

where /3 is the emissivity index, Tjj is the dust temper- 
ature in Kelvin, So is an amplitude factor, and S(v) is 
in mJy. In principle, it would be better to do the analy- 
sis without making an assumption about the opacity of 
the dust. However, the quality of the 350/im data is not 
sufficient to dist inguish r from (1 — e ~ r ) in fitting SEDs 
to the data (see IKovacs et all [200 6c). Furthermore, the 
number of data points is generally insufficient to fit more 
than two parameters; assuming the dust is optically thin 
eliminates the source size from the fit. We return to this 
issue in Sec. [6l 

The luminosity of each source was calculated by inte- 
grating the best fit SED: 

L = (2.5 x HT 11 ) 4ttD 2 l J S u dv (2) 

where Dl is the luminosity distance in Mpc and L 
is the luminosity in L Q . In Table 4, the Lpin, val- 
ues are for integrals of the SEDs over rest frequencies 
that correspond to observer frame wavelengths of 42.5- 
122.5 /im, using the definition of the FIR band given by 
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iSanders and Mirabel I (|1996f ); the Ldust values are for in- 
tegrals over all frequencies. Since the shortest wavelength 
data we consider is A reS f = 58 /im, we assume that all 
the emission is coming from dust grains. In a few cases, 
there may be contributions from other emission mecha- 
nisms to the longest wavelength data (see but these 
points generally do not drive the fit. The dust temper- 
ature, Td, and (3 are highly correlated and cannot be 
separately determined with t he few data points t ypically 
available for fitting SEDs (see lBeelen et al. II2006D . In fit- 
ting SEDs to the observed data points, we held (3 fixed, 
with the value f3 — 1.5. This value is approximately that 
determined in da ta fits where there w as sufficient infor- 
mation to vary (3 (jBeelen et al. 1 12006) . The value chosen 
does not have a dramatic effect on the calculated lumi- 
nosity, which is determined largely by the observed 350 
/im point. In our data fits, changing the fixed value of (3 
by ± 20% changed the luminosity by about +20/ — 50%. 
Overall, we regard our calculated luminosities typically 
to be accurate to within a factor of 2-3. 
Dust masses given in Table 4 were calculated from 

M d = (4.8 x 10 14 ) • (3) 

(1 + z) K{v rest )B(v rest ,T D ) 

where Md is in M Q , B v is the Planck function, and 
K>u = kq{v i 'IT 'H z) 13 is the dust mass absorption coeffi- 
cient. For Td we used the results of the SED fit to the 
data. Similarly to calculated luminosities, the calculated 
values of Md depend on the value of (3 that is assumed. 
For v rest we have chosen u rest — 0.35(1 + z) THz, the rest 
frequency for an observed wavelength of 850 /im, where 
the spectrum is solely dust emission. We set k = 1.3 
m 2 kg" 1 at 1 THz by int e rpola ting in the table given by 
Ossenkop f and Henning I (|1994f) . 

Table 4 also lists the star formation rate (SFR) cal- 
culated from L F IR using SFR(M e yr" 1 ) = (1.8 x 
10- 10 )L FIR (L Q ) (|Kennicutt I I199H. L dust /M dust , and a 

minimum radius, R^in™ j f° r a source of the dust radi- 
ation, assumed to be a disk seen face-on. We calculated 
Rmin m °y assuming the source to be optically thick at 
Vrest = 0.86(1 + z) THz, the rest frequency for an ob- 
served wavelength of 350^m. While some of the energy 
to heat the dust in the qu asars may be supplied by accre- 
tion onto the black hole (jGranato et al.|[1996D . we have 
attributed all of it to star formation. Consequently, the 
star formation rates we derive could be overestimated. 

4. DISCUSSION OF INDIVIDUAL SOURCES 

4.1. Detections 

LB QS 18 is an optically identified quasar from the 
survey of iFoltz et ait (11989)) detected i n CO(3-2) with 
the IRAM Interferometer (jlzaakj USSl ■ LBQS0018 is 
radio quiet. 

SMM J02396 w as detected in CO(3-2) emission by 
iGreve et all (|2005l) using the IRA M Interferometer. The 
source is a ring galaxy identified bv lSoucail et al. 1 (|1999l ) 
from HST images. They measured the redshift of z = 
1.062 and determined the mag nification factor b y the 
cluster Abell 370 to be n = 2.5. iSmail et all (|2002D refer 
to the source as J02399-0134, rather than J02396-0134 
as it has come to be labeled. 

RX J0911.4 is a ROSAT source identified as a mini- 



broad absorption line quasar bv lBade et al. I (0.997). The 
quasar is strongly l ensed (// = 22) by a galaxy at z 
= 0.8 (jKneib et al.ll200l . RX0911.4 was detected in 
CO (3-2) at the Ow ens Valley Radio Observatory by 
Hai nline et all (|2004h . The quasar is radio quiet. 

SMM J14011 was the second SCUBA source to be 
detected in a molecular emission line, na mely CO (3- 
2) at the Owens Valley Radio Observatory (|Fraver et al.l 
1999). It has been intensively stud ied since then in CO 
lines (see lDownes fc Solomonll200l : the CI(l-O) line has 
been detected (jWeifi et al. 1120051 )7 The source is lensed, 
but the lens magnification is very uncertain and model 
depen dent, ranging from /i = 5 — 25 (|Downes &: Solomon! 
[2001 . 

4.2. Tentative Detections 

SMM J02399 i s a hyperluminous infrared galaxy 
(jlvison et al.l 119981 ). the first SCUB A source to be de - 
tected in molecular line emission, bv lFraver etHl (|1999l ) 
at the Owens Valley Radio Observatory in CO(3-2). 
From higher angular resolution ob s ervat ions with the 
IRAM Interferometer, iGenzel et al.l (|2003l ) argued for a 
model with a very large (r > 2 kpc) disk, a possibility 
that awaits still higher angular resolution confirmation. 
Our detection at 350 /im of 29 ± 9 mJy is somewhat 
incon sistent with the high value at 450 /im (jlvison et alJ 
119981 ). (Fig. 3), even considering the large errors in each. 

SPS S 0338 is a high-r edshift (z = 5.03) quasar discov- 
ered bv lFan et all ([1999). The CO (5-4) line was detecte d 
with the IRAM Interferometer bv lMaiolino et alJ (|2007|) . 
The radio continuum has not been detected . This source 
was pr eviously detected at 350 ^m by iWang et al. I 
(|2008al) . with S 350lim = 17.7 ± 4.4 mJy. Our value of 
29 ± 10 mJy is consistent within the uncertainties. 

MG0751 is a quasar detected in CO(4-3) emission by 
iBarvainis et alJ (|2002al) using the IRAM Interferometer. 
The quasar is highly lense d, with magnification of 16.6 
(jBarvainis fc IvisCTl2002bT ). This is our weakest detec- 
tion, with a S/N ratio of only 2.3, weaker than the 450 
/jm detection. 

SMM J 09431 is a submillimeter ga laxy detected in 
CO (4-3) emission bv lNeri et all ((20031) using the IRAM 
Interferometer. 

5. TESTING THE IDEA OF BASIC UNITS 

While high redshift starbursts are quantitatively more 
extreme than anything in the local universe, some of 
their intensive properties, such as the infrared luminos- 
ity per HCN line luminosity, are similar to those of lo- 
cal starbursts an d even cluster-forming den s e clumps in 
our o wn Galaxy (|Gao and Solomon~ll2004aH bl; IWu et al~l 
2005). This similarity has led to the suggestion that we 
can understand them in terms of a large number of basic 
units of star formation, with these units being patterned 
on the well-studi ed massive dense clumps in our Galaxy 
(|Wu et al. 1120051 ) . This model is similar in spirit, though 
not in detail, to a model suggested by ICombes et al.l 

We use the data on this sample of EMGs to make a 
reality check on this proposal. The idea is simple. Start- 
ing with some mean properties of massive, dense clumps 
in our Galaxy, we divide the luminosity of the EMG into 
Ndump units and consider the consequences. 
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We take as our s ample of Galact ic massive, dense 
clumps the survey of IWu etldl (|2009f) . which brings to- 
gether data on 5 molecular lines and infrared emission 
from over 50 regions, including some of the most lumi- 
nous regions in our Galaxy, such as W31, W43, W49, 
and W51. It is important to include these, as the distri- 
bution of luminosities is strongly skewed to lo wer values, 
and w e seek a mean value. To be consistent, IWu et all 
( 2009) calculated the luminosity of the massive clumps 
with the same method as that used for high-z galaxies. 
The mean luminosity of this sample is (4.7 ± 1.2) x 10 5 
Lq . IWu et aLl (|2005l ) noted that the ratio of infrared to 



line luminosity of these clumps matches that of starburst 
galaxies only for clumps with a luminosity above about 
10 4 ' 5 Lq. If we restrict the average to those clumps, the 
mean value is (6.3±1.6) x 10 5 L Q . For simplicity, we will 
take a value of 5 x 10 5 L© as the standard value. Table 5 
shows the number of such average clumps that would be 
needed, calculated simply from N c i ump = Lpm/5 x 10 5 . 
As explained in we adopted a value of 1 x 10 5 L Q 
for SMM J02396, with the corresponding increase in the 
number of clumps. Values oiN c i ump range from 8.8 x 10 5 
to 2.8 x 10 7 . 

Next, we ask if this large number of units would rea- 
sonably fit into the space available. If not, the basic unit 
model could not work, and we would have to turn to 
models with still denser dust, perhaps in a torus around a 
black hole. We seek a minimum size estimate. The most 
economical packing is spherical, but this packing would 
likely result in high optical depth, even at submillime- 
ter wavelengths, so we also consider a two dimensional 
packing, in which the units are placed in a disk, only 
one unit thick. The truth should be somewhere between 
these two estimates. The minimum radii are then given 

1 /3 

by Rsphere = r ciumpN c [ ump for the spherical packing and 



Rdl 



r du mp N^ mj) f or th e disk packing. Based on 

the sample in IWu et all (|2009h . (r c i ump ) = 1.13 ±1.09 pc 
as measured by the size of the HCN 1-0 emission. For 
simplicity, we take a standard value of 1.0 pc. Sizes for 
spherical packing range from 95 to 300 pc, while those 
for disk packing range from 940 to 5300 pc. These sizes 
are not unreasonable, though the concept of a 5 kpc disk 
filled uniformly with massive dense clumps may strain 
credulity. We emphasize that these are minimum size 
estimates and that close packing of such clumps is very 
artificial. This calculation is only a reality check and not 
a serious model. We present such a model in the next 
section. 



6. MODELING THE UNITS 

The usual approach to modeling continuum emission 
from high-z starbursts is to fit an optically thin, sin- 
gle temperature, modified black body. The modification 
is to assume an opacity that varies with wavelength as 
K v = ^(Are/XA/Are/) • This is the procedure we have 
followed also, as described in Section [5] For observers of 
Galactic star formation, this seems a quaintly anachro- 
nistic procedure. Even a single region has a distribution 
of temperatures, and opacities can be approximated by 
a power law only at long wavelengths. In galaxies with 
large redshifts, the rest wavelengths observed in the sub- 
millimeter are in the far-infrared, where dense clumps 



can be optically thick. These effects are commonly mod- 
eled by radiative transfer calculations that use realistic 
opacities, place a luminosity source in the center of a 
clump, and correctly account for optical depth effects. 

By taking the basic unit idea a step further, we can ex- 
plore the effects of these complications by running radia- 
tive transport models for an entity representing the mean 
clump. The observational properties that constrain the 
models are, as usual, the flux densities. To perform ra- 
diative transfer in the rest frame of the emitting galaxy, 
we scale all wavelengths to the rest frame and we scale 
the flux densities according to 



S v (clump) — 



5 x 10 5 



£fir \ 10k P c / M 



Surest), (4) 



where (x is the lens magnification factor. 

The choice of 10 kpc for the standard distance to the 
clump from the observer is arbitrary and chosen just for 
convenience in comparison to clumps in our Galaxy. 

Because we keep the luminosity and size fixed at the 
values adopted above (L = 5 x 10 5 L Q and r c i ump = 
1.0 pc), the only variables are those that describe the 
mass and distribution of matter in the clump. Models of 
continuum radiation from massive, dense clumps show 
that that they are well modeled by power law density 
distri butions (e.g., iMueller et al. I [20021 : iBeuther et al. I 
l200l : 

n(r) =n f (r)(r/r f )- p (5) 

For example, IMueller et al. I (|2002l ) fit radiative transfer 
mode ls to data for 31 sourc es and found (p) = 1.8 ± 0.4, 
while IBeuther et al. I (pOOl) found (p) = 1.6 ± 0.5 in the 
inner regions, steepening in the outer regions. We can 
explore the effects of different choices for p. We adopt 
single power laws over the whole clump. For connection 
to previous work, we take 77 = 1000 AU. The mass of 
the clump is then determined by the combination of n/ 
and p. The optical depth through the clump depends on 
the density distribution and the assumed inner radius. In 
modeling Galactic clumps, the longest wavelength point 
is generally used to constain the mass because it is the 
most likely to be optically thin. Thus, for a given value of 
p, rif is adjusted until the predicted flux density matches 
the observations. For quasars, there could be contamina- 
tion by free-free or synchrotron emission at the longest 
wavelength. The only modeled source with a peculiar 
data point at long wavelengths is MG0751 and we ignore 
that point in our fit. 

We use the radiative transfer code by lEgan et al. I 
(1988) as modified by us to produce outputs convenient 
to our p urpose. We use the opac i ties ta bulated in col- 
umn 5 of lOsscnko pf and Henningl ([1994) , commonly re- 
ferred to as OH5. The OH5 opacities have been shown 
to reproduce well many observations of massive dense 
clumps. There is some evidence that they result in un- 
derestimates of the mass by a factor of about 2, but this 
is still uncertain and undoubtedly varies from region to 
region. 

The luminosity in the models is represented by a single 
star at the center of the mass distribution. The lumi- 
nosity available for heating the dust is set to the stan- 
dard value of 5 x 10 5 Lq. Depending on the effective 
temperature of the star, its actual luminosity must be 
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greater to account for photons used to ionize an HII re- 
gion. With no information to constrain the nature of the 
forming star cluster, we simply set the effective tempera- 
ture to 44,000 K, roughly equivalent to an 05.5 star, for 
which half the luminosity is in Lyman continuum pho- 
tons. This choice is largely irrelevant, as many studies 
have shown that the stellar photons below the Lyman 
limit are rapidly absorbed and degraded to longer wave- 
lcngt hs, so that the initial s pectrum is rapidly erased 
(e.g.. Ivan der Tak et al.lll999f l . The models also include 
heating from the outside by an interstellar radiation field 
(ISRF). The ISRF is based on studies of the Milky Way 
near us, as shown in lEvans et al. I ([2001). Since the ra- 
diation field is likely to be much higher in these galaxies, 
we multiply the standard ISRF by a factor of 100 as our 
default value. We will explore the effect of changing this 
value. The models also require an inner radius to the 
envelope. We adjust this value so that the temperature 
there is in the range of 1000 to 2000 K, because dust will 
evaporate at higher temperatures. The choice of inner 
radius has no impact on the mass of the envelope or the 
flux density at long wavelengths, but it has a major ef- 
fect on the optical depth. The optical depth can strongly 
affect the SED at shorter wavelengths, but the effect at 
the wavelengths we are modeling is minor. The models 
result in values for M gas , the clump mass, assuming a gas 
to dust ratio of 100 by mass, the ratio of the luminosity 
to mass, and the optical depth at all wavelengths, which 
we characterize by tioo, the optical depth at 100 /im. 

6.1. Example: Application to RXJ0911.4 

As an example and to explore the effects of parame- 
ter choices, we consider the case of RXJ0911.4, a highly 
lensed quasar at z — 2.8. It has data available from 350 
/zm to 3000 /im, or from 92 to 790 /im in the rest frame. 

Figure[4]shows the SED with the data and several mod- 
els for this source, which are defined in Table 5. Figure 
[S] shows the run of density and dust temperature in each 
of the models. In Model 1, we assume p = 2, and n/ 
was adjusted to match the longest wavelength point. It 
also matches the SHARC-II data, which corresponds to 
a rest-frame wavelength of 92 /im. However, it over- 
estimates all the data at intermediate wavelengths. If 
the data at rest-frame 790 /im has a contribution from 
sources other than dust emission, the mass could be over- 
estimated. There was also a lower value obtained at 
the same (observed) wavelength of 3000 /im (§4}, which 
might suggest variability. The rest of the models have a 
lower value of n/ and produce less emission at long wave- 
lengths. They provide a better compromise between the 
various observations. In Model 3, we change p to 1.6, and 
Model 4 has p = 1.0. We adjust n/ in each model to keep 
the predicted flux density about the same at rest-frame 
790 /im, about 60% of the (higher) observed value. As is 
apparent, the value of p has little effect on the SED at 
the wavelengths with observations. However, the value 
of p does affect the optical depth. The values of optical 
depth at 100 /im in the rest frame are given in Table 5. 
These show that the usual assumption of optically thin 
emission is highly dubious if the inner radius is as small 
as our standard value. Data at shorter wavelengths are 
needed to constrain the optical depth, though the as- 
sumption of spherical symmetry then becomes an issue. 
For lower values of p, higher masses are required to pro- 



duce the same flux density at 790 /im because the lower 
values of optical depth result in more of the mass being 
at lower temperatures. However, the difference in mass 
is about a factor of two, as shown in Table 5. 

We also show Model 5, which is the same as Model 
2 except that we increase the ISRF to 5000 times the 
local value. This produces an upturn in the tempera- 
ture at large radii, where the ISRF takes over, and the 
SED is affected at the shorter wavelengths. In particu- 
lar, the model now is much too strong at the rest-frame 
wavelength of the SHARC-II observations. The total lu- 
minosity from this model, based on integrating the SED, 
more than doubles. Since this extra luminosity is coming 
from outside and would represent luminosity from adja- 
cent star forming regions or exposed luminous stars, it is 
just a redistribution of the total luminosity of the galaxy. 
However, Figure |4] shows that a sufficiently well sampled 
SED could begin to distinguish these situations. 

For the other sources, we restrict modeling to ISRF of 
100 and p — 2. We use the variation in Ld U st/M gas found 
for the various values of p for RXJ0911.4 to estimate an 
uncertainty of 40% on derived values of L dust/ M gas- 

6.2. Two Extreme Cases 

Having established the basic dependencies on parame- 
ters for RXJ0911.4, we discuss only a few salient features 
of models for the other sources. 

SDSS0338, a quasar at z = 5.0, is the most extreme 
source in our sample. With a Lfir = 14 x 10 12 L Q and 
a SFR — 2500 M Q yr _1 , it exemplifies the most lumi- 
nous sources known. As noted in Table 4, it has the 
warmest SED, with To = 56 K. The modeling of this 
source, shown in Figure [6l is limited by having only 3 
data points, none longer than 200 /im in the rest frame. 
Data at longer wavelengths would be the most valuable 
addition to the data base on this source. Table 5 shows 
that the high value of To really implies a high luminosity 
to mass ratio, about five times that of RXJ0911.4. With 
the usual interpretation of this ratio as a "star formation 
efficiency" , a more useful way to look at the shape of the 
SED is in terms of Ld ust /M gas than in terms of a single 
temperature. 

SMM J02396, at z = 1.06, is the coldest (T D = 22 K) 
and least luminous object in our sample, with Lfir, — 
2 x 10 11 L Q and a modest star formation rate, SFR = 37 
M Q yr _1 , only about 10 times that of the Milky Way. In- 
terestingly, it cannot be fitted with the standard proce- 
dure. As shown in Figure[3 the inferred flux densities for 
a standard clump are so high that models fail to fit. To 
get even close to the longest wavelength point, extremely 
high fiducial densities are needed, which cause huge opti- 
cal depths. The resulting cold SED indeed tries to repro- 
duce the low value for Td from the optically thin fitting, 
but it falls far below the shorter wavelength data. Essen- 
tially, a source characterized by such a cold SED cannot 
have a luminosity as high as 5 x 10 5 L . The lower panel 
of Figure [7] shows the result of lowering the luminosity 
of the standard clump to 1 x 10 5 L Q . Now the flux den- 
sities are also five times lower (see equation 4) and can 
be fitted with a reasonable model. The resulting ratio of 
luminosity to mass is very low. In this galaxy, the mean 
clump would be less luminous and less star-forming than 
the mean clump in our galaxy. The galaxy is detected at 
350 /im only because it has a huge amount of dense gas 
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compared to the Milky Way. 

This galaxy is a good candidate for the kind of cir- 
rus m odel proposed by lEfstathiou &: Rowan-Robinson! 
(2003), where the dust responsible for the far-infrared 
emission is far from the starburst. The fact that it is 
a ring galaxy (|Soucail et al. I [1999) is consistent with 
this kind of picture, but current observations cannot tell 
whether the far-infrared emission is extended on the scale 
of the 4'.' 8 diameter ring. 

7. DISCUSSION 
7.1. Relations 

The value of L dust /M gas can be related to a depletion 
time for the dense gas and is often described as a star 
formation "efficiency" . Values from the radiative trans- 
fer models range from 1.9 to 283, with a mean of 105 and 
a median of 78. A similar calculation for dense clumps 
in the Milky Way, as defined by HCN maps, yields sim- 
ilar values: in a range of 3.3 to 398, the mean value is 
117 ± 20, with a median of 97. This similarity indicates 
that modeling extreme starbursts as an ensemble of mas- 
sive clumps similar to those in the Milky Way leads to a 
consistent result. 

In the absence of detailed models, how should one in- 
terpret the dust temperature derived from fitting the 
SED? There is a strong correlation between the fitted 
dust temperature and the luminosity to mass ratio for the 
mean clump derived from the radiative transport models 
(Figure HJ . A least-squares fit in log space to the points 
with uncertainties in both axes yields the following: 

log(L dust /M gas ) = (-4.55 ± 1.54) + (4.11 ±0.97)log(T D ) 

(6) 

We have left SMM J02396 out of the fit, as it may repre- 
sent a galaxy dominated by cirrus emission. It is plotted 
in Figure [8] and clearly deviates from the fit to the other 
sources. The value of 4.11 ± 0.97 can be understood in 
terms of the Stefan-Boltzmann law, though we caution 
that Td is, at best, an average over a large range of real 
temperatures. A similar relatio n is found for m assive, 
dense clumps in the Mil ky Way dWu et aljl2009h . 

Using the relation from fcennicutt I (|1998f ) between star 
formation rate and far-infrared luminosity, we can com- 
pute a depletion time for the dense gas from 

t dep (y) = 5.6 x 10 9 {L dust /M gas )- x « 2.0 x lO 14 !^ 411 

(7) 

which ranges from 3.0 x 10 9 yr for SMM J02396 down to 
2.0 x 10 7 yr for SDSS0338. The mean value is 4.7 x 10 s 
yr, or 1.1 x 10 8 yr, excluding the value for SMM J02396. 

The bottom plot in Figure[8]shows the log of the optical 
depth at 100 /im from a single core versus the value of 
dust temperature from the simple fit to the SED. There 
is a strong anti-correlation. (We have again left SMM 
J02396 out of the fit, but it lies on the relation in this 
case.) Galaxies with SEDs suggesting warmer dust can 
be modeled with less opaque clumps as more luminosity 
is derived from a smaller amount of mass. 

7.2. Comparison of Dust and CO 

The dust masses we have calculated from the simple, 
optically thin isothermal fit (Table 4) assume an opacity 
given by models of dust in dense clouds in the Milky Way 



(|Ossenkopf and Henning |["l994f) . The gas masses emerg- 
ing from the more detailed models and given in Table 5 
make the same assumption about opacity and also as- 
sume a gas to dust ratio of 100, as is usually assumed in 
the Milky Way. Both these assumptions could in princi- 
ple be quite wrong for high-z EMGs. Masses computed 
from CO luminosities depend on a conversion to molecu- 
lar hydrogen based on models of local ULIRGS, which is 
about a fact or of 6 times lower than for gala xies like the 
Milky Way (|Solomon fc Vanden Bout1[2005l ). Thus sub- 
stantial uncertainties attend these mass estimates and 
consistency checks are worthwhile. 

We have compared the masses derived from the simple 
fit to those from the detailed radiative transport models 
by using a gas to dust ratio of 100 to convert the former 
into gas masses. The mean ratio of simple estimate to 
model mass is 0.64 with surprisingly small dispersion. 
Thus, the simple fit gives a good estimate of the mass, 
even when the assumptions are invalid. As long as there 
is a data point at wavelengths long enough to be optically 
thin and in the Ray leigh- Jeans limit, a simple fit provides 
an estimate that can agree with those from more detailed 
models. If the models are correct, the masses from simple 
fits should be scaled up by a mere 1.56. 

Takin g the gas mass estimates b a sed on CO luminos- 
ity from iSolomon fc Vanden Bout I (|2005f) . we find that 
the mean ratio of the mass from CO to that from the 
models is 0.60, with a bit more dispersion (minimum is 
0.15 for SMM J02396 and maximum is 1.50 for SMM 
J14011). The agreement between these two estimates is 
even more surprising, as there are many reasons that they 
could differ. The conversion from CO luminosity to mass 
is uncertain and depends on inte rnal cloud conditions 
like density and tem perature (e.g.. iDickman et al.lll986l 
ISolomon et al.l[l987t ). The dust opacity model could be 
different or the gas to dust ratio could differ. For exam- 
ple, the two mass estimates could be reconciled if the gas 
to dust ratio is 167, which is probably wit hin the r a nge o f 
uncertainty, even for the Milky Way (see iDrainej (2003) 
for discussion of abundance issues in grain models) . The 
general consistency despite all these possible differences 
suggests that neither the standard conversion from CO 
luminosity to gas mass nor the gas to dust ratio is likely 
to be far off. The less satisfying possibility, that both 
conversions are wrong and the errors cancel out, cannot 
be ruled out. 

Considering the possible differences in radiation envi- 
ronment, metallicity, etc., the consistency between the 
masses from CO and from dust emission is quite remark- 
able. We have used the submillimeter opacity values 
appropriate for coagulated, icy grains in dense molec- 
ular cores. The submillimeter opacity of the dust grains 
that characterize the diffuse interstellar medium in the 
Milky Way is less by a factor of 4.8, which would result 
in masses from dust emission that are, on average, larger 
than those inferred from CO by a factor of 3. Assuming 
that the CO conversion factor is exactly correct, the dust 
properties in EMGs would appear to lie between those 
of the diffuse ISM in the Milky Way and those of dense 
cores, but closer to the latter. 

The remarkable thing about this comparison is that 
the masses from the models include only quite dense gas 
(typically n > 10 4 cm" 3 , see Fig. 4), while the CO 
emission could easily arise in gas of much lower density. 
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The agreement suggests that essentially all the molecular 
gas in these EMGs is in dense clumps. 

7.3. Comparison to Other Models 

iNarayanan et al.1 (|2009af ) have recently proposed a 
model for the formation of high redshift submillimeter 
galaxies as the result of a starbursts triggered by ma- 
jor mergers in massive halos. They modeled the time 
evolution of such mergers and found that the most lumi- 
nous observed submillimeter galaxies could be modeled 
as 1:1 mergers in the most massive dark halos. As the 
two galaxies spiral in toward the merger, SFRs of 100- 
200 M© yr -1 can be triggered. During final coalescence, 
they predicted a brief burst (5 x 10 yr) of very high 
SFR (~ 2000 Mq yr- 1 ). These predictions are broadly 
consistent with our numbers for SFR (Table 4) and our 
mean depletion time of (1.1-4.7) x 10 8 yr. To make the 
most luminous submillimeter galaxies, they need to as- 
sume that all stellar clusters with ages less than 10 Myr 
are embedded in the material that is forming them, con- 
si stent with the assumptio ns of our modeling. 

INarayanan et al.1 (|2009b| ) have discussed the interpre- 
tation of CO lines from such galaxies. They warn that 
the mass estimates from higher-J CO lines may under- 
estimate the gas mass by a factor up to 2. This would 
be consistent with our result above, but we believe that 
the sources of uncertainty discussed in §7.2 are equally 
important. 

It is also interesting to compare our results 
to predictions of quite diff e rent kinds of models. 
lEfstathiou fc R owan-R obinsonl (|2003f ) modeled two of 
our sources with their model of cirrus dominated far- 
infrared emission. In SMM J02399, our data are consis- 
tent with their model, but in SMM J14011, their model 
prediction is about a factor of 10 too low. This result is 
generally consistent with the To values in Table 5: SMM 
J02399 has T D — 29 K, within the limit of 30 K that the 
dust achieves in their models; SMM J14011 has T D = 42 
K. On this basis, about half the sources in our sample 
could be consistent, with SMM J02396 as a particularly 
good case. Of course, models in which some of the radi- 
ation arises in cirrus and some in dust surrounding the 
forming stars could be relevant. Since the cirrus model 
is the logical extreme of model 5 of RXJ0911.4, where 
we increased the ISRF, further constraints on the SED 
could help to constrain such models. 

While our data do not bear directly on models using 
radiation from an AGN to heat a dusty torus, the general 
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consistency of our models with dust emission from star 
formation bolsters the case for the embedded starburst 
models. 

8. SUMMARY 

We have observed 15 galaxies with redshifts from 1 
to 5 at 350 /im. Four have detections at levels above 
5(7, while four have detections of lower significance. Far- 
infrared luminosities range from 2 x 10 11 to 1.4 x 10 13 
Lq, and inferred star formation rates range from 37 to 
2500 Mq yr" 1 . From fits to optically thin, isothermal 
emission with an opacity index, (3 — 1.5, characteristic 
dust temperatures range from 22 to 56 K. 

Because the dust is unlikely to be optically thin and 
isothermal, we consider a picture in which the star form- 
ing regions are composed of a large number (N c i ump = 
0.9 — 30 x 10 6 ) of dense clumps, each with a luminosity 
equal to 5 x 10 5 L Q , roughly the mean value for massive 
star forming clumps in the Milky Way. A crude calcu- 
lation of the minimum size needed to pack such a large 
number of clumps into a galaxy does not rule out such 
models. 

Radiative transport models of standard clumps are 
then used to match the observed SEDs, converted into 
the rest frame and scaled to a single dense clump at a 
distance of 10 kpc, for ease of comparison to Milky Way 
clumps. These models show that the individual clumps 
are likely to be quite opaque at far-infrared wavelengths, 
but that the simple fits do capture the total mass of emit- 
ting dust quite well. The differences between sources lie 
primarily in the ratio of luminosity to mass, which is 
commonly taken as a measure of star formation "effi- 
ciency" in extragalactic studies. Values of 2 to 283 are 
found for Ld U st/M gas . This ratio shows a strong correla- 
tion with the value of To from the isothermal, optically 
thin fit, resulting in Ldust/M gas cx To 4 ' 11 . The depletion 
timescales for dense gas range from 3 x 10 9 y for SMM 
J02396 down to 2.0 x 10 7 y for SDSS0338. 
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TABLE 1 
Source List and Observations 



Source 


R.A. 
(J2000) 


Decl. 
(J2000) 


Redshift 
(z) 


Observation 
Date 


Zenith 
Angle 


T225GHz 
(neper) 


Integration 
(hr) 
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nn-21 -27 3n 
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2.62 
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20 


0.06 


1.0 










12/06 


30 


0.05 


2.0 










10/07 


25 


0.06 


2.0 


SMM J02396 


02-39-56 60 


-01 -34-26 6 


1.06 


12/06 


30 


0.04 


2.0 










10/07 


25 


0.05 


1.0 


SMM J02399 


02:39:51.90 


-01:35:58.8 


2.81 
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25 


0.05 
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snssnrns 

O J— I kj kj \J vJ v> O 


03-38-29 31 
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50 
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Gao, Y., & Solomon, P.M. 2004a, ApJ, 606, 271. 
Gao, Y., & Solomon, P.M. 2004b, ApJS, 152, 63. 
Genzel, R., Baker, A. J., Tacconi, et al. 2003, ApJ, 584, 633. 
Granato, G. L., Danese, L., & Franceschini, A. 1996, ApJ, 460, 
Lll 

Greve, T. R., et al. 2005, MNRAS, 359, 1165. 

Hainline, L. J., Scoville, N. Z., Yun, M. S., Hawkins, D. W., 

Frayer, D. T., & Isaak, K. G. 2004, ApJ, 609, 61 
Ivison, R. J., Smail, I., Le Borgne, et al. 1998, MNRAS, 298, 583. 
Ivison, R.J., Smail, I., Barger, et al. 2000, MNRAS, 315, 209. 
Izaak, K.G. 2004, private communication. 
Kennicutt, R. C, Jr. 1998, ARA&A, 36, 189. 
Kneib, J.-P., Co hen, J. G., fe Hjor th, J. 2000, ApJ, 544, L35. 
Kovacs, A. 2006, arXiv:0805.3928v2. 
Kovacs, A. 2006b, PhD thesis, Caltech. 

Kovacs, A., Chapman, S. C, Dowell, et al. 2006c, ApJ, 650, 592. 
Leong,M.M. 2005, URSI Conf. Sec, J3-J10,426 
Maiolino, R., et al. 2007, A&A, 472, L33. 

Mueller, K.E., Shirley, Y.L., Evans, N.J., II, et al. 2002, ApJS, 
143, 469. 

Narayanan, D., Hayward, C. C, Cox, T. J., Hern quist, L., 

Jonsson, P., Younger, J. D., & Groves, B. 2009, arXiv:0904.0004 

Narayanan, D., Cox, T. J., Hayward, C, Younger, J. D., & 
Hernquist, L. 2009, arXiv:0905.2184 

Neri, R., et al. 2003, ApJ, 597, L113. 

Ossenkopf, V. & Henning, T. 1994, A&A, 291, 943. 

Priddey, R.S., Isaak, K.G., McMahon, R. G., et al. 2003a, 
MNRAS, 339, 1183. 



Priddey, R. S., Isaak, K. G., McMahon, et al. 2003, MNRAS, 344, 
L74. 

Rowan-Robinson, M. 2000, MNRAS, 316, 885 
Sanders, D.B., & Mirabel, I.F. 1996, ARA&A, 34, 749. 
Smail, I., Ivison, R.J., Blain, A. W., et al. 2002, MNRAS, 331, 
495. 

Solomon, P.M. & Vanden Bout, P.A. 2005, ARA&A, 43, 677. 
Solomon, P. M., Rivolo, A. R., Barrett, J., & Yahil, A. 1987, ApJ, 
319, 730 

Soucail, G., Kneib, J.P., Bezecourt, J., et al. 1999, A&A, 343, 
L70. 

Spergel, D.N., et al. 2007, ApJS, 170, 377. 

van der Tak, F. F. S., van Dishoeck, E. F., Evans, N. J., II, et al. 

1999, ApJ, 522, 991. 
Wang, R., Carilli, C.L., Wagg, J., et al. 2008a, AJ, 135, 1201. 
Wang, R., Carilli, C.L., Wagg, J., et al. 2008b, ApJ, 687, 848. 
Weifi, A., Henkel, C, Downes, D., et al. 2003, A&A, 409, L41. 
Weifi, A., Downes, D., Henkel, C, et al. 2005, A&A, 429, L25. 
Weifi, A., Downes, D., Neri, R., Walter, F., Henkel, C, Wilner, 

D. J., Wagg, J., & Wiklind, T. 2007, A&A, 
Wu, J., Evans, N.J. II, Gao, Y., et al. 2005, ApJ, 635, L173. 
Wu, J, Dunham, M. M.; Evans, N. J., II, Bourke, T. L., Young, 

C. H., 2007, AJ, 133, 1560 
Wu, J., Evans, N. J., II, Shirley, Y. L., & Knez, C, in prep. 



TABLE 2 

Source Detections 



Source 


Flux Density 


Noise 


S/N 




(mJy) 


(mJy) 




LBQS0018 


32 


5 


6.4 


SMM J02396 


51 


6 


8.5 


SMM J02399 


29 


9 


3.1 


SDSS0338 


29 


9.5 


3.1 


MG0751 


36 


16 


2.3 


RXJ0911.4 


150 


21 


7.1 


SMM J09431 


22 


6.6 


3.3 


SMM J14011 


75 


10 


7.5 



TABLE 3 

Measured Flux Densities of Detected Sources 



Source 


350/im a 


450fim 


750fj,m 


850^tm 


1200/im 


1300/wn 


1350^tm 


3000/im 


Ref b 




(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 




LBQS0018 


32±5 






17.2±2.9 










1 


SMM J02396 


51±6 


42±10 




11±1.9 










2 


SMM J02399 


29±9 


69±15 


28±5 


26±3 






5.7±1.0 




3 


SDSS0338 


29±10 






11.9±2.0 


3.7±0.3 








1.1 


MG0751 


36±16 


71±15 




25.8±1.3 






6.7±1.3 


4.1±0.5 


5.6 


RXJ0911.4 


150±21 


65±19 




26.7±1.4 




10.2±1.8 




1.7±0.3 


6 


SMM J09431 


22±7 






10.5±1.8 




2.3±0.4 






7,8 


SMM J14011 


75±10 


41.9±6.9 




14.6±1.8 






2.5±0.8 




9,10 



a Al l the 350^m fluxes are f rom this work. 

b 1. IPriddev et al.l lt2003al1. 2. ISmail et al.l ||2002T). 3. Ilvison et all (Il998h. 4. ICarilli et afl f2WS). 5. IBarvaims et~aTI 
H2002al). 6. IBarvainis fc Ivisonl l|2002bl ). 7. ICowie et al.l (|2002l ). 8. INeri et alj H2003l )79. Ilvison et al.l H20001) . 10. 
IDownes k, S olomon (20031). 



TABLE 4 

lNTRiNSic a Derived Parameters: Dust Temperatures, Luminosities, Masses; Star Formation Rates 



Source 


Lens 


Tdust 


Lfir 


Ldust 


850m 

JlU s t 


SFR 


L -dust 1 'M dust 


mm 




Mag. 


(K) 


(1O 12 L ) 


(10 12 L ) 


(10 s M ) 


(Moyr- 1 ) 


(10 4 L /M ) 


(pc) 


LBQS0018 


1 


28 


2.9 


4.5 


9.6 


520 


0.53 


1710 


SMM J02396 


2.5 


22 


0.2 


0.1 


1.9 


37 


0.21 


300 


SMM J02399 


2.5 


29 


1.8 


2.6 


5.8 


320 


0.45 


1030 


SDSS0338 


1 


56 


11 


23 


2.6 


2500 


8.9 


775 


MG0751 


17 


31 


0.4 


0.6 


0.9 


72 


0.75 


470 


RXJ0911.4 


22 


37 


0.8 


1.0 


0.4 


140 


2.5 


435 


SMM J09431 


1.2 


39 


4.3 


5.5 


3.1 


770 


1.8 


830 


SMM J14011 


5-25 


42 


2.2-0.4 


2.9-0.6 


0.7-0.15 


400-90 


4.1-40 


445-200 



a Corrected for the effects of lensing. 
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TABLE 5 
Model Parameters 



Source 


N clump 


R(sphere) 


R(disk) 


Model 


n f 


V 


noo 


Mgas 


Ldust/Mgas 






(pc) 


(PC) 




(cm-' 3 ) 






(Mq) 


(Lq/Mq) 


t Rrmnm s 


fi Q v 1 f)6 

o.y x iu 


loU 




A 
1 


l.o X 1U 


9 n 


I 




94 


SMM J02396 a 


2.0 x 10 6 


125 


1414 


5 


3.3 x 10 9 


2.0 


125 


5.3 x 10 4 


1.9 


SMM J02399 


3.6 x 10 6 


153 


1910 


2 


1.0 x 10 9 


2.0 


60 


1.6 x 10 4 


31 


SDSS0338 


2.8 x 10 7 


302 


5300 


1 


1.1 x 10 s 


2.0 


11 


1.8 x 10 3 


283 


MG0751 


8.8 x 10 5 


95 


940 


3 


1.0 x 10 9 


2.0 


60 


1.6 x 10 4 


30 


RXJ0911.4 


1.6 x 10 6 


116 


1260 


1 


6.7 x 10 s 


2.0 


39 


1.1 x 10 4 


47 










2 d 


4.0 x 10 s 


2.0 


23 


6.4 x 10 3 


78 










3 


9.2 x 10 7 


1.6 


11 


8.8 x 10 3 


57 










4 


6.6 x 10 6 


1.0 


2 


1.1 x 10 4 


47 










5 


4.0 x 10 s 


2.0 


23 


6.4 x 10 3 


113 b 


SMM J09431 


8.5 x 10 6 


203 


2920 


2 


3.0 x 10 s 


2.0 


18 


4.8 x 10 3 


104 


SMM J14011 c 


4.4 x 10 6 


162 


2090 


2 


1.5 x 10 s 


2.0 


9 


2.4 x 10 3 


208 



a This model uses a clump luminosity of 1 X 10 5 Lq, five times less thant the standard value, used for all other galaxies. 

Model 5 is the same as Model 2, except that the ISRF has been multiplied by 5000 instead of 100; the extra luminosity is 
all from the external radiation. 
c Assumes a lens magnification of 5. 



This model was used in correlations and statistics. 
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a (J2000) „ (J2000) 




a (J2000) 



Fig. 1. — Images of 350 lira emission for four galaxies with detections at 5a or better. The black circle at the lower right show the size 
of beam. The colors of images from blue to red indicate increasing intensity for the pixels. The intensity contours start from the 2a level, 
with a step of 2a. The la levels for LBQS0018, SMMJ02396, RXJ0911.4 and SMMJ14011 are 5, 6, 21, 10 mjy within a 20" aperture, 
respectively. The black cross in each panel marks the position of the CO peak. The green circles indicate the region within which detection 
is reliable. Emission towards the edge is not reliable because of decreased integration time and should be ignored. 



12 




04 s 02 
a (J2000) 



a (J2000) 



Fig. 2. — Images of 350 fim emission for four galaxies with detections at less than 5a. The black circles at the lower right show the 
size of beam. The colors of images from blue to red indicate increasing intensity for the pixels. The intensity contours start from 1.5<x, 
with a step of 1.5a, except for MG0751, in which the contour levels are 1.5a and 2.2a. The la levels for SMMJ02399, SDSS0038+0021, 
SMMJ09431 and MG0751 are 9, 9.5, 9, 16 mjy within a 20" aperture, respectively. The black cross in each panel marks the position of 
the CO peak. The green circles indicate the region within which detection is reliable. Emission towards the edge is not reliable because of 
decreased integration time and should be ignored. 
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Fig. 3. — The SEDs of all eight sources in the observed frame. 
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Fig. 4. — The SED of the mean clump in RXJ0911.4 in the rest frame frequencies and flux densities at a distance of 10 kpc. The models 
are described in the text and in Table 4. 
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Fig. 5. — The density and dust temperature are plotted versus radius for the models of the mean clump in RXJ0911.4. The models are 
described in the text and in Table 4. 
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Fig. 6. — The SEDs for the standard clump, placed at 10 kpc, in six galaxies with detections are plotted as points with errorbars. The 
flux densities are shown at the rest wavelengths and the flux densities are scaled using equation 4. The solid lines are the best-fitting 
models from the radiative transfer modeling, as indicated in Table 4. 
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Fig. 7.— (Top) The SED for the standard clump (L = 5 X 10 5 L Q ), placed at 10 kpc, in SMM J02396. No model was able to fit the flux 
densities. (Bottom) The SED for a less luminous standard clump (L = 1 X 10 5 L©), placed at 10 kpc, in SMM J02396. Now a model can 
provide a reasonable fit. 
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1 1.2 1.4 1.6 1.8 2 
log T d (K) 

Fig. 8. — (Top) The log of the luminosity to mass ratio from the radiative transport models is plotted versus the log of the dust temperature 
derived from the isothermal, optically thin fit. The solid line is the least squares fit to the data given by equation 6. The point at very 
low Tp is SMM J02396, which has been excluded from the fit. (Bottom) The log of the optical depth at 100 /im from the best model is 
plotted versus the log of the dust temperature derived from the isothermal, optically thin fit. The solid line is the least squares fit to the 
data. The uncertainties assigned to the optical depth were 100%, to roughly indicate the dependence on model parameters. The point at 
very low To is SMM J02396, which has been excluded from the fit. 



